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The Monoceros radio loop: temperature, brightness, spectral index and 
distance 
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In this paper we estimated the temperatures and brightnesses of the Monoceros radio loop at 1420, 820 and 408 MHz. 
Linear spectrum is estimated for mean temperatures versus frequency between 1420, 820 and 408 MHz. The spectral 
index of Monoceros loop is also obtained. The brightness temperatures and surface brightnesses of the loop are computed 
using data taken from radio-continuum surveys at the three frequencies. The spectral index of the loop is also obtained 
from T - T plots between 1420 - 820, 1420 - 408 and 820 - 408 MHz. The obtained results confirm non-thermal origin 
of the Monoceros radio loop. 
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1 Introduction 

The Monoceros filamentary loop nebula was suggested to 
be a supernova remnant (SNR) by Davies 1963 on the basis 
of 237 MHz observations. This shell source of radio emis- 
sion was suspected already by Davies 1963 to be a SNR, 
while Gebel & Shore 1972 developed this idea into a more 
detailed study. It was considered as an object similar to ma- 
jor loops when Spoelstra 1973 included it in his study of 
galactic loops as supernova remnants expanding in the local 
galactic magnetic field. 

Graham et al. 1982 found that the constellation of Mo- 
noceros is remarkably rich in extended Galactic radio sour- 
ces. Large parts of it have been mapped in the radio con- 
tinuum over a wide range of frequencies. Monoceros Neb- 
ula can be found in a catalog of Galactic SNRs listed as 
G205.5H-0.5 (Green 2004, Green 2006). It is S-type remnant 
(shell remnant), characterized by diffuse, shell-like emis- 
sion with steep radio spectrum and size of '--^220'. Shell type 
supernova remnants are believed to be particle accelerators 
to energy up to a few hundred TeV, and it is shown that 
Monoceros loop is a good candidate for acceleration of par- 
ticles (Fiasson et al. 2007). 

The aim of this paper is to calculate the average bright- 
ness temperatures and surface brightnesses of the Mono- 
ceros radio loop at 1420, 820 and 408 MHz and to study 
how these results compare with previous result (Urosevic 
& Milogradov-Turin 1998; Milogradov-Turin & Urosevic 
1996) and with current theories of supernova remnant evo- 
lution. These theories predict that SNRs are at radio fre- 
quencies predominantly non-thermal sources which are sp- 
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reading inside of the hot and low density bubbles made by 
former supernova explosions or by strong stellar winds (see 
Salter 1983; McKee & Ostriker 1977 and references there- 
in). 

Spectrum (temperature versus frequency) have been pl- 
otted and this result is used to determine spectral index of 
the Monoceros loop. 



2 Analysis 



2.1 Data 

In this paper we used observations from several radio-conti- 
nuum surveys given in "Flexible Image Transport System" 
(FITS) format which are available on MPIfR's Survey Sa- 
mpler ("Max-PIanck-Institut fiir Radioastronomie", Bonn). 
This is an online service (http://www.mpifr-bonn.mpg.de/| 
survey.html), which allows users to pick a region of the 
sky and obtain images and observed data (in FITS for- 
mat) at a number of wavelengths. User can choose a coordi- 
nate system, projection type and a survey. The radio con- 
tinuum surveys at 1420, 820 and 408 MHz provided the 
database for computing brightness temperatures (7b)- The 
1420-MHz Stockert survey (Reich & Reich 1986) has the 
angular resolution of 35', the 820-MHz Dwingeloo survey 
(Berkhuijsen 1972) 1°.2 and the 408-MHz all-sky survey 
(Haslam et al. 1982) 0°.85. The corresponding observations 
are given at the following rates for both I and b: ^ at 1420 
MHz, ^ at 820 MHz and ^ at 408 MHz. The effective 
sensitivities are about 50 mK T^, 0.20 K and about 1 .0 K, 
respectively. 
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2.2 Method 

We extracted observed brightness temperatures from FITS 
format into ASCII data files and after that, these data files 
have been processed by our software, i.e. we have developed 
several programs in C and FORTRAN in order to obtain 
results presented in this paper. 

The area of Monoceros loop is very difficult to deter- 
mine precisely due to great influence of background radia- 
tion and superposed external sources, such as Rosette Neb- 
ula (Graham et al. 1982; Aharonian et al. 2004, Aharonian 
et al. 2007; Kim et al. 2007). Therefore, some authors made 
very rough estimates for the shape of this loop, represent- 
ing it by fitted circles (see e.g. Aharonian et al. 2004; Kim 
et al. 2007). In order to make better estimates for the loop 
boundaries, we analyzed temperature profiles (see Figure[Tli 
for different values of Galactic latitude (b) between -6° and 
6°, and for Galactic longitude (0 from 212° to 198°. As 
one can see in Figure [T] there are two temperature peaks. 
The lower one corresponds to the total brightness temper- 
ature of the loop (including background radiation), while 
the higher one corresponds to superposed external sources. 
We denoted the minimum and maximum brightness tem- 
perature of the lower peak by Tmin and T^ax, respectively. 
These temperature limits enable us to distinguish the loop 
from background and also from external sources, and their 
values are given in Table [T] Moreover, using the tempera- 
ture profiles we found that the Monoceros loop is located in 
this sky region: I € [210°, 200°] and b G [-6°, 5°]. 

The area of Monoceros loop is enclosed with bright- 
ness temperature contours (see Figure |2]). These contour 
lines correspond to the minimum and maximum brightness 
temperatures which define its borders (see Table [Til, and 9 
contours in between. We have subtracted the background 
radiation, as well as the superposed radiation, in order to 
derive the mean brightness temperature of the SNR alone 
(see Borka 2007 and Borka, Milogradov-Turin & Urosevic 
2008). 

For deriving temperatures over the Monoceros loop, the 
areas used for the individual spurs were obtained from the 
radio continuum maps at three frequencies: 1420 MHz (Re- 
ich & Reich 1986), 820 MHz (Berkhuijsen 1972) and 408 
MHz (Haslam et al. 1982). The areas over which an average 
brightness temperature is determined at each of the three 
frequencies are taken to be as similar as possible within the 
limits of measurement accuracy. Temperature limits Tj^in 
and Tmax, given in Table[Tl are the lower and the upper tem- 
perature values: Tj^in is the lower temperature limit between 
the background and the spur, while Tmax is the upper tem- 
perature limit between the spur and superposed confusing 
sources. In this manner, background radiation was consid- 
ered as radiation that would exist if there were no spurs. If 
the value of Tmin (Tmax) is changed by some small value 
(approximately 0. 1 K for 1420 MHz and more for other fre- 
quencies), the brightness contours become significantly dif- 
ferent. For evaluating brightness temperature of the back- 
ground, we used all measured values below Tmin, inside the 



corresponding intervals of I and b, and lying on the outer 
side of a loop. The value of Tb is approximately constant 
near a loop. For evaluating the brightness temperature of a 
loop including the background, we used all measured values 
between T^nin and Tmax inside the corresponding regions of 
I and b. Mean brightness temperature for loop is found by 
subtracting the mean value of background brightness tem- 
perature from the mean value of the brightness temperature 
over the area of the loop. After deriving the temperature, we 
derived the surface brightness by: 

S = (2A:z/Vc')Tb (1) 

where k is Boltzmann constant and c is the speed of light. 
Results are given in Table [T] 

The surface brightnesses of SNRs must be above the 
sensitivity limit of the observations and must be clearly dis- 
tinguishable from the Galactic background emission (Green 
1991). Therefore, the data from the fainter parts of the loops 
(which are very low surface brightness SNRs) are not taken 
into account because it is very difficult to resolve them from 
the background. On the other hand, this would significantly 
reduce brightness of entire loop and there is a general trend 
that fainter SNRs tend to be larger (Green 2005). For evalu- 
ation brightness temperatures over the spurs we had to take 
into account background radiation (Webster 1974). Borders 
enclosing the spurs are defined to separate the spur and its 
background. 

The measured data have different resolutions for differ- 
ent frequencies (see §2.1), and therefore in order to obtain 
r — T plots the data are retabulated so the higher resolution 
maps are convolved to the resolution of the lowest resolu- 
tion map. In that way we convolved data at 1420 and 408 
MHz to 0°.5 X G°.5 resolution, which is the sampling rate of 
the 820 MHz survey. These retabulated data are presented in 
Figure|4]for the following frequencies: 1420 MHz (top left), 
820 MHz (top right), 408 MHz (bottom left). Then, for each 
frequency pair we used only the common points (with the 
same {I, b)) which belong to the loop area at both frequen- 
cies (see Figure IH bottom right). In that way we reduced 
loop area to the same area for different frequencies. The ob- 
tained T — T plots for three pairs of frequencies (between 
1420-820, 1420-408 and 820-408 MHz) enabled calculating 
the spectral indices. For each of the three frequency pairs, 
by interchanging the dependent and independent variables 
we have obtained two j3 values for each pair and the mean 
value of these fit results is adopted as the radio spectral in- 
dex, as suggested in Uyaniker et al. 2004. 

3 Results and discussion 

The average brightness temperatures and surface brightnes- 
ses of Monoceros radio loop at the three frequencies are 
given in Table [T] 

Spectrum shows how temperature depends on the fre- 
quency: 

log T = log if-/? log t/ (2) 
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Fig. 1 Temperature profiles at 1420 MHz for Galactic longitude from 212° to 198° and for the following values of 
Galactic latitude: l.°5 (top left), 1° (top right), 0° (middle), -2° (bottom left) and -3° (bottom right). Temperatures are 
given in K, and Galactic longitudes in degrees. 
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Fig. 2 The map of a region in Monoceros, in new galactic coordinates (1, b). Contours of the brightness temperatures Tb 
are plotted. Top left: the 1420 MHz map with contours of Tb given in units of K. The contours are plotted every 0.04 K 
starting from the lowest temperature of 3.8 K up to 4.2 K. Top right: the 820 MHz map with Tb in K. The contours are 
plotted every 0.19 K starting from 8.8 K up to 10.7 K. Bottom left: the 408 MHz map with Tb in K. The contours are 
plotted every 1.1 K starting from 36 K up to 47 K. Bottom right: temperature colorbars for 1420, 820 and 408 MHz, all 
in K. 
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Table 1 Temperatures and brightnesses of Monoceros radio loop at 1420, 820 and 408 MHz. 



Frequency 


Temperature limits 


Temperature 


Brightness 


(MHz) 




(K) 


(10-22 W/(m^ Hz Sr)) 


1420 


3.8, 4.2 


0.18 ±0.05 


1.09 ±0.30 


820 


8.8, 10.7 


0.90 ± 0.20 


1.85 ±0.40 


408 


36, 47 


5.2 ± 1.0 


2.63 ± 0.50 



where if is a constant, and /3 is radio spectral index. 

The spectrum was generated using mean temperatures at 
three different frequencies. This best-fit straight line spec- 
trum enables calculation of spectral index as negative value 
of the line's direction coefficient. The spectrum is shown in 
Figure [3] Relative errors of the measurements AlogT — 
rp'f^-j^Q are presented by error bars, where AT are the corre- 
sponding absolute errors given in Table[T] Obtained value f3 
= 2.70 ± 0. 14 (greater than 2.2) confirm non-thermal origin 
of Monoceros loop emission. The value for the brightness 
temperature spectral index of the Monoceros loop is rather 
steep (about 2.7). This is at the high end of the spectral in- 
dex distribution for SNRs as suggested in Clark & Caswell 
1976. 

From the relation {T — T plot): 
r.i/T,2 = iiyi/iy2y (3) 
P can be determined as 

/? = log(ai2)/log(j^i/j^2), (4) 

where ai2 is direction coefficient of line Ti,2{T,yi). 

Because of the influence of superposed sources to the 
loop and because it is hard to estimate precise loop's bor- 
ders, we present how spectral indices vary for different b. 
Variations of spectral indices between the three frequency 
pairs obtained from T-T plots, distributed over Galactic lat- 
itudes, are given in Tables |2]|4] For some Galactic latitudes, 
there is a large dispersion of the points in T — T graphs so 
they are not suitable for calculation of spectral index; we 
did not calculate j3 for them and it is denoted with The 
averaged values of radio spectral indices are given in Table 
121 Examples of T - T plots between 1420-820, 1420-408 
and 820-408 MHz for longitude 6 = 1 ° are given in Figure|5] 
The average value of spectral index from T — T is < /3tt > 
= 2.63 ± 0.30. It can be noticed that this value agrees well 
with the corresponding value obtained from spectrum, as 
expected (see Uyaniker et al. 2004). Then we calculated 
mean value of spectral index between 1420, 820 and 408 
MHz (regarding spectrum and T — T graphs) < (3 > = 2.66 
± 0.20. With this < /? > we reduced to 1000 MHz (Ta- 
ble |6l) according to relation: 

SlGHz/S.GHz- Wl)^''"'^. (5) 

In Graham et al. 1982, the value (3 = 2.47 ± 0.06 of radio 
spectral index of this loop was derived from the radio con- 
tinuum spectrum between two frequencies: 111 and 2700 
MHz, and we calculated spectrum at three frequencies and 
also derived T — T graphs which represent different method. 



Table 2 Variations of spectral indices between 1420 and 
820 MHz obtained from T-T plots, distributed over Galactic 
latitudes. 



latitude (°) 


/3(1420 - 820) 


/3(820 - 1420) 


-5.0 






-4.5 






-4.0 






-3.5 






-3.0 


1.92 ±0.82 


4.54 ± 0.82 


-2.5 


3.86 ± 0.46 


4.84 ± 0.46 


-2.0 


3.47 ±0.31 


4.07 ±0.31 


-1.5 


2.96 ± 0.73 


5.02 ± 0.73 


-1.0 


1.96 ±0.55 


3.31 ±0.55 


-0.5 


2.61 ±0.39 


3.36 ± 0.39 





2.39 ± 0.22 


2.73 ± 0.22 


0.5 


2.49 ± 0.28 


3.09 ± 0.28 


1.0 


2.42 ± 0.26 


2.85 ± 0.26 


1.5 


2.74 ± 0.27 


3.13 ±0.27 


2.0 


2.62 ± 0.26 


3.02 ± 0.26 


2.5 


1.59 ± 0.41 


2.69 ± 0.41 


3.0 


1.26 ±0.43 


241 ± 0.43 



Distances to the SNRs can be inferred from positional 
coincidences with HI, HII regions and molecular clouds, 
OB associations, or pulsars or from measuring optical ve- 
locities and proper motions. Where there is no direct dis- 
tance determination, estimates can be made for shell rem- 
nants by utilizing the radio surface brightness-to-diameter 
relationship (S — D) (Case & Bhattacharya 1998 and refer- 
ences therein). The mean surface brightness at a specific ra- 
dio frequency, T,^, is a distance-independent parameter and, 
to a first approximation, is an intrinsic property of the SNR 
(Shklovskii 1960). 

Using a set of 172 SNRs as calibrators, for which re- 
liable distance values were determined, Urosevic 2002 has 
constructed the following relation between the surface brig- 
htness (S) and diameter (D): 

EiGHz = 2.76 X lO-^^D-^ °^ (6) 

This relation, combining the Galactic and extragalactic 
SNRs (master relation), also includes four main Galactic ra- 
dio loops (Loops I - IV). 

It is shown (Odegard 1986) that there is diffuse emission 
with the filaments located in the western half of Monoceros 
loop (overlapping with Rosette Nebula) and in the eastern 
part (dust cloud L1631). There is also a peak in the 2700 
MHz emission in the eastern edge for which Graham et al. 
1982 suggested to be a HII region. In Aharonian et al. 2007 



www.an-journal.org 



© WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



6 



V. Borka Jo\ano\ ic and D. Uroscvic: The Monoccros radio loop 




log V 

Fig. 3 Monoceros loop spectram: temperature versus frequency, for three measurements - at 408, 820 and 1420 MHz. 



Table 3 Variations of spectral indices between 1420 and 
408 MHz obtained from T-T plots, distributed over Galactic 
latitudes. 



latitude (°) 


/3(1420 - 408) 


/3(408 - 1420) 


-5.0 


1.92 ± 0.24 


2.32 ± 0.24 


-4.5 


2.14 ±0.33 


3.15 ±0.33 


-4.0 


2.35 ±0.12 


2.62 ± 0.12 


-3.5 






-3.0 


2.48 ±0.19 


3.03 ±0.19 


-2.5 


2.91 ± 0.13 


3.14 ±0.13 


-2.0 


2.58 ± 0.11 


2.78 ±0.11 


-1.5 


2.55 ±0.14 


2.82 ± 0.14 


-1.0 


2.55 ± 0.08 


2.65 ± 0.08 


-0.5 


2.69 ± 0.07 


2.76 ± 0.07 





2.58 ±0.10 


2.75 ±0.10 


0.5 


2.39 ±0.16 


2.78 ±0.16 


1.0 


2.40 ±0.19 


2.82 ± 0.19 


1.5 


2.16 ±0.22 


2.65 ± 0.22 


2.0 


2.19 ±0.25 


2.80 ± 0.25 


2.5 




3.23 ± 0.89 


3.0 







Table 4 Variations of spectral indices between 820 and 
408 MHz obtained from T-T plots, distributed over Galactic 
latitudes. 



latitude (°) 


/3(820 - 408) 


/3(408 - 820) 


-5.0 


2.38 ± 0.51 


3.29 ±0.51 


-4.5 


1.39 ± 0.89 


4.20 ± 0.89 


-4.0 


1.71 ±0.49 


3.22 ± 0.49 


-3.5 


1.12 ±0.56 


2.74 ± 0.56 


-3.0 


1.87 ± 0.33 


2.74 ± 0.33 


-2.5 


1.57 ± 0.28 


2.07 ± 0.28 


-2.0 


1.47 ± 0.25 


1.95 ±0.25 


-1.5 


1.24 ± 0.35 


2.07 ± 0.35 


-1.0 


2.23 ±0.31 


2.87 ± 0.31 


-0.5 


2.44 ±0.17 


2.67 ± 0.17 





2.49 ± 0.10 


2.58 ±0.10 


0.5 


2.23 ± 0.12 


2.36 ±0.12 


1.0 


2.22 ±0.18 


2.50 ±0.18 


1.5 


1.91 ± 0.25 


2.36 ± 0.25 


2.0 


1.96 ±0.18 


2.29 ±0.18 


2.5 




3.21 ± 0.72 


3.0 







it is said that the complex Monoceros loop/Rosette Neb- 
ula region contains several potential sources of very-high- 
energy 7 ray emission and that the interaction of the SNR 
with a compact molecular cloud is possible. Because of pos- 
sible influence with molecular clouds and HII region, we 
also take S — Z) relation from Arbutina et al. 2004 adjusted 
for molecular clouds: 



SiGHz = 1.1 X IQ-'^D 



-3.5 



(7) 



They have taken 14 shell- type Galactic SNRs associated 
with large molecular clouds. 

Using the relations derived by Urosevic 2002 and Ar- 
butina et al. 2004, we computed diameter of the loop and 
then we computed distance as: 



r = D/{2smQ), 



(8) 



Table 5 The average radio spectral indices between 1420, 
820 and 408 MHz obtained from T-T plots. 



Frequency pair (MHz) 



1420 - 820 
1420 - 408 
820 - 408 



2.98 ± 0.40 

2.62 ± 0.20 
2.29 ± 0.34 



Table 6 Brightnesses of Monoceros radio loop reduced to 
1000 MHz using /? = 2.66 ± 0.20. 



Frequency Brightness at 1000 MHz 
(MHz) (10-2^ W/(m^ Hz Sr)) 



1420 
820 
408 



1.37 ±0.50 
1.63 ±0.30 
1.45 ±0.03 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.org 




Astron. Nachr. / AN () 




1 



I 



I 



I 



I 



208 



206 204 
1 



202 200 



□□□□□□□□□□□□□□□□□□ 
□□□□□□□□□□□□□□ □□□□ 

□□□□□□□□□□□□a □ 
□□□□□□□ □□□□ □□ 



□□□□□□ □ 
□□□□□□□□□ 

□□□□□□□□ 

□□□□□ 
□□□□□ 
□□□□□ 
□□□□□ 
□□ □ 



□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□□ 
□□□□ □□□□□□ 



□□□□□□ □□□□□□□□□□□□ 
□□□□□□□□□□□□ □ □□ □□ 

□□□□□□□□□□□□ □□□□□□ 
□□□□□□□□□□□ □ □□□□□ 
□□□□ □□□□□ 



— I — ' — I — ' — I — ' — I — ' — I — ' — I — 

210 208 206 204 202 200 



1 



Fig. 4 The data retabulated to 0° .5 x 0°.5 resolution, for the following frequencies: 1420 MHz (top left), 820 MHz (top 
right) and 408 MHz (bottom left). Bottom right: the common pixels which belong to the loop area at all frequencies. 



with angular radius (9) taken from Urosevic & Milogradov- 
Turin 1998. The results are given in Table|2l 



able T, - D relation by Arbutina et al. 2004. If not, then it 
is more suitable to use relation by Urosevic 2002. 



It can be noticed that there is the influence of molecular 
cloud on Monoceros loop (Graham et al. 1982, Aharonian 
et al. 2007). If that influence constantly existed during the 
loop's evolution, then for distance calculation it is more suit- 



For the distance to this loop, mean optical velocity sug- 
gests 0.8 kpc, and low frequency radio absorption suggests 
1.6 kpc (Graham et al. 1982,Odegard 1986, Green 2006). 
With our calculated radio spectral index (3 = 2.66 ± 0.20, 
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Fig. 5 Top left: T-T plot between 1420 and 820 MHz for Monoceros loop at latitude b = 1°, in longitude range [210°, 
200°]. Top right: the same as top left, but between 1420 and 408 MHz. Bottom: the same as top left, but between 820 and 
408 MHz. 



Table 7 Diameters D (pc) and distances r (pc) of Mono- 
ceros radio loop derived from the S — Z? relation given by 
Urosevic 2002 and by Arbutina et al. 2004. 



Relation 


D(pc) 


r (pc) 




119 ± 18 


1630 ± 250 




92 ± 14 


1250 ± 190 



we calculated distance to this loop using our derived bright- 
nesses. The results are given in Table|7] 

4 Conclusions 

In this paper we calculated the brightness temperatures and 
surface brightnesses of the Monoceros radio loop at 1420, 
820 and 408 MHz. Linear spectrum is estimated for mean 
temperatures versus frequency between 1420, 820 and 408 
MHz. It is the first time that the brightness temperatures of 



Monoceros loop are calculated at 820 and 408 MHz fre- 
quencies from the observational data. We sampled much 
more points (more than 1 000) at 1420 MHz than in previ- 
ous papers (95 points) (Urosevic & Milogradov-Turin 1998, 
Milogradov-Turin & Urosevic 1996). Also, the brightness 
temperature is now derived using the different method. The 
temperature of this radio loop at 1420 MHz is in good agree- 
ment with the result obtained in Urosevic & Milogradov- 
Turin 1998. 

We present the radio continuum spectrum of the Mono- 
ceros loop using average brightness temperatures at three 
different frequencies. As it can be seen from Figure[3] given 
linear fit provides reliable spectral index. Also, we present 
the T — T plots which enables also calculation of spectral 
index. 

The effective sensitivity of the brightness temperatures 
are: 50 mK for 1420 MHz, 0.2 K for 820 MHz, and about 
1.0 K for 408 MHz. The most precise measurements (the 
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least relative errors) are in case of 1420 MHz, so positions of 
the brightness temperature contours of the loop are the most 
realistic for this frequency. The brightnesses of Monoceros 
loop at 1420 and 820 and 408 MHz are in good agreement 
when reduced to 1 GHz. 

With our derived brightnesses, we calculated new di- 
ameters and distances to this loop at the three frequencies: 
1420, 820 and 408 MHz and then estimated some average 
distance. We used empirical Y, — D relations for supernova 
remnants by Urosevic 2002 and Arbutina et al. 2004. The 
estimated distance of the Monoceros radio loop is in good 
agreement with the earlier results (Guseinov et al. 2004; 
Green 2006: 1.6 kpc). 

The spectral index analysis confirms non-thermal origin 
of Monoceros radio loop. 
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